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Abstract 

The time variation of fundamental mass scales can have profound cosmo- 
logical implications. We investigate a particular model of crossover quintessence 
which is compatible with all present cosmological observations. This model 
can also reconcile the reported time variation of the fine structure constant 
from quasar absorption lines with the bounds from archeo-nuclear physics and 
tests of the equivalence principle. 
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Dynamical mass scales 



The recent confirmation of dark energy by the high precision KMAP data [1] 
challenges for an answer to the question if dark energy is constituted by a static 
cosmological constant or by dynamical quintessence [2]. A dynamical dark energy 
may shed light on the profound question of the role and origin of mass scales in 
physics. It may also lead to a cosmological time variation of dimensionless funda- 
mental parameters like the electromagnetic fine structure constant. 

In quantum field theory or statistical physics the dimensionless couplings depend 
on length scales [3]. This effect is perhaps most striking in QCD where the running 
gauge coupling generates the mass scale for the hadrons. The modern view of the 
renormalization group [4] associates this running with an "integrating out" of fiuc- 
tuations as one looks at the effective physics on larger and larger length scales. The 
fiuctuation effects can be treated in a conceptual simple and elegant form in the con- 
text of an exact functional renormalization group equation for the cff'cctive average 
action [5]. It describes for both renormalizable and nonrenormalizable theories the 
dependence of the fiuctuation effects on an effective infrared cutoff. The physical 
length scale which cuts off the fiuctuations in the infrared may correspond to the 
expectation value of a scalar field [6] and therefore become dynamical. We explore 
here the hypothesis that the time evolution of scales persists until the present epoch 
in the universe [2] . This would have striking effects for cosmology, the time variation 
of fundamental "constants" and tests of the equivalence principle. 

For the renormalization of the standard model of particle physics one usually 
assumes implicitely the existence of an ultraviolet scale - typically the Planck scale 
or a grand unification scale. Similarly, for solids the ultraviolet scale is given by 
the lattice distance or for gases and liquids by the molecular size. In contrast, we 
explore here the concept of a unified theory where the action does not involve any 
explicit fundamental mass scale ^. Even without knowing the fundamental unified 
theory the experience with quantum physics lets us expect two types of possible 
mass scales in such a scenario. First, some intrinsic mass scales may be generated 
by the running of dimensionless couplings, similar to the hadron masses in QCD. 
Second, the expectation values of scalar fields induce dynamical mass scales. This 
second type of scales can change in the course of the cosmological evolution. 

We explore here the scenario that the Planck mass or Newtons constant is dy- 
namical. In contrast to the early formulation of this idea by Dirac and Jordan [7], 
[8] we postulate [9], [2] that the particle masses are dynamical as well, such that 
in a first approximation the ratio between the nucleon mass and the Planck mass, 
rrin/Mp, remains fixed, and similar for the weak scale M^/Mp, the electron mass 
■rrie/Mp etc. In this first approximation all mass scales are proportional to a scalar 
field X, i.e. Mp — x, m„ = /i„x or me = heX- This changes profoundly both the 

^In string theories the apparent mass scale corresponding to the string tension can be replaced 
by the expectation value of the dilaton, thereby making this scale dynamical as well. 
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conceptual status and the observational consequences as compared to scenarios with 
fixed particle masses [8, 10]. 

In fact, in absence of intrinsic mass scales generated by the running dimension- 
less couplings the effective action (after integrating out all quantum fluctuations) 
would exhibit an exact dilatation symmetry, corresponding to a multiplicative rescal- 
ing X ^ (^X- If '^^ neglect in a first approximation the intrinsic mass scales, any 
nonzero value of the scalar field x would break the exact global dilatation symme- 
try. In consequence of the spontaneous breaking of scale symmetry most particles 
would acquire a mass ~ x- Another consequence would be a Goldstone boson ~ In^ 
which has only derivative interactions - the dilaton. In this approximation of an ex- 
act spontaneously broken dilatation symmetry no interesting consequences for late 
cosmology are expected [9, 2]. However, the impact of the dilaton on cosmology 
changes dramatically if the effect of the possible intrinsic mass scales is taken into 
account. In terms of the dilatation symmetry the intrinsic mass scales refiect an 
"anomaly" induced by the quantum fluctuations. The scalar fleld x can now con- 
stitute a dynamical dark energy or quintessence [2] which is relevant for the present 
universe. Because of its important cosmological consequences the "pseudo-dilaton" 
In X has been named the cosmon. 

We present here a particular version of crossover quintessence [11] with two 
intrinsic mass scales generated by running dimensionless couplings. The flrst mass 
scale Xc indicates a crossover behavior in the derivative interactions of the pseudo 
Goldstone boson, without introducing a potential. More concretely, it affects the 
kinetic term of the cosmon and indirectly the location of an ultraviolet fixed point 
for the running grand unified gauge coupling. The second scale m is analogous to the 
hadron scale Aqcd in QCD. It is an infrared scale where the running grand unified 
gauge coupling grows large. This scale m induces a mass term for the cosmon field 
~ w?x^. Similar to the ratio Kqcd/ Mp the ratio m/xc can be exponentially small. 
We will see that for \n{m/xc) ~ 1/138 a reahstic cosmology arises, with quintessence 
responsible for about 70% of the present energy density of the universe. Our model 
turns out to be compatible with all present cosmological tests. 

Furthermore, this model predicts that the fundamental "constants" depend on 
time [9, 2]. In fact, the running of a typical dimensionless coupling like the fine struc- 
ture constant ctem can be written in the functional form Q;em(?^/x^, X^/a*^) where // 
is an arbitrary renormalization scale. For fixed x the dependence on the squared 
momentum is governed by the usual /9-function of the standard model. On the 
other hand, for fixed q^/x^ the variation of x is governed by the running of the 
grand unified gauge coupling ax- The cosmological time dependence of x results 
in a time dependence of the fine structure costant! We show that for reasonable 
/^-functions for the running grand unified gauge coupling our model can reconcile 
the claimed time variation from the observation of quasar absorption lines [12] with 
severe bounds from archeo-nuclear physics [13, 14] and tests of the equivalence prin- 
ciple [15]. 
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Crossover Quintessence 



Crossover quintessence (CQ) can be characterized by a recent increase of the 
fraction in homogeneous dark energy = Ph/Pcr from a small but not negligible 
value for early cosmology to a value Qh ~ 0.7 today. A small early value (say 
flh = 0.01) typically results from a cosmic attractor solution [2, 16, 17, 18] ("tracker 
solution" ) independently of the detailed initial conditions. For such models of "early 
quintessence" the dark energy decreases in early cosmology at the same pace as ra- 
diation. Dark energy and radiation have therefore always been of a comparable 
magnitude and excessive "fine tuning" can be avoided [2, 19]. The amount of early 
quintessence is accessible to observation by nucleosynthesis [2] [20], the cosmic mi- 
crowave background (CMB) [21] or structure formation [22, 23]. The recent crossover 
to a domination of the universe by dark energy (Q/j ^ 0.7) may be related to prop- 
erties of the underlying model ^ or triggered by matter domination [24] or structure 
formation [25]. On the level of the equation of state CQ corresponds to a crossover 
from a value Wh ~ 1/3 during the radiation dominated universe to a substantially 
negative value for z < 0.5. 

In this note we elaborate on a recent proposal [11] that CQ is connected to the 
existence of a conformal fixed point in a fundamental theory. This will serve as a 
natural starting point for a discussion of the variation of couplings. We assume that 
for low momenta (g^ -C x^) the yet unknown unified theory results in an effective 
grand unified model in four dimensions. There the coupling of the cosmon field x 
to gravity and the gauge fields is described by an effective action 

S = Jd':r^!^-^x'R+l{S{x)-Q)d'xd,X 

+m V + ^F'^^F.J ■ (1) 



The running of the effective grand unified gauge coupling gj^ = / Zp is determined 
by the x-dependence of Zp and we assume a fixed point behavior for ax = fi'x/47r 

^""^ -hax-hal-h^a\S{5). (2) 



dhix 

We take here a small value of [fte] such that the last term only leads to a small shift 
in the fixed point value 

- ^ - ^- (3) 

(We choose positive 62 and hi and S{5) will be specified below.) For ax{x) > Oix,* 
the grand unified gauge coupling increases with decreasing x- Similar to QCD it 
grows large at some nonperturbative scale ("confinement scale") that we associate 

^This typically requires a modest "tuning" of parameters at the level somewhat below 1%. 
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with m. Our model therefore has an "infrared scale" m which is generated by 
dimensional transmutation from the running of the gauge coupling. We assume 
that m determines the mass of the cosmon. Neglecting the effects of the dilatation 
anomaly reflected by m the cosmon potential would vanish. 

For X ^ the gauge coupling runs towards its fixed point (3) as x increases. 
At the fixed point the effective action would have an effective dilatation symmetry if 
the running of 5{x) could be neglected. We will assume, however, that the running 
of S is instable towards large x (£■ > 0) and postulate 

Ps{S) = E5^ (4) 



dlnx 

For small enough 5^ = 5(x = m) the dimensionless coupling 5 stays small for a large 
range m < % < Xc- For x ^ ?7t. we may neglect the small mass m for the evolution 
of the dimensionless couplings. The evolution in this range is then governed by the 
vicinity to a conformal fixed point [11] at 5 = 5* = , ax = ax,*- In the same 
approximation our model corresponds to a flat direction ^ in the effective potential, 
consistent with conformal symmetry. All cosmon interactions are therefore given 
by its derivative coupling to the graviton and gauge fields. We note that the fixed 
point behavior bears a certain resemblance with the "runaway dilaton" . [26] 

Nevertheless, as x increases further an effective "ultraviolet scale" Xc character- 
izes the region where S grows large according to the solution 

The ratio between the ultraviolet and infrared scales turns out exponentially huge 

At the crossover scale Xc the flow of the couplings witnesses a crossover from the 
vicinity of the conformal flxed point (small 6) to another regime (perhaps a different 
fixed point) for large 5. We will see that the crossover behavior in the present epoch 
of the cosmological evolution precisely corresponds to the crossover from small to 
large S for x in the vicinity of Xc- Realistic cosmology obtains for E6i ~ 1/138. 

The precise behavior of f3s in the region of large 6 will only be of secondary 
relevance for past and present cosmology provided 6 grows large enough near Xc- 
For the present note we investigate first a scenario where the increase of S with x is 
unbounded while 5{x) is defined for all x, namely 

*In string theories such flat directions are often called "moduli". In presence of additional 
matter fields, e. g. scalar and fermions, the flat direction in the space of scalar fields corresponds 
to constant ratios of all masses [11]. 
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(We take = 0.05 and observe that cosmology shows no strong dependence on the 
value of Js. This holds provided that Js is small enough such that the difference 
between eqs. (4) and (7) matters only in the region of large 5. At the end of this 
note we will compare this scenario with an alternative where (3s exhibits a fixed 
point for large 5.) For very large 5 in the region x ^ Xc the solution obeys now 
^(x)/^(Xc) ^ ix/Xc)^^'^^ ■ We observe that for large 5 a standard scalar kinetic term 
obtains for a rescaled scalar field % ~ X^^^^(x) ~ In this language the 

coupling to gravity vanishes for x — > oo according to y^R ~ 5~^x^R. Our first 
scenario for the behavior at large 5 can therefore also be viewed as a fixed point in 
the cosmon-gravity coupling at 5^^ = 0. 

In the sense of a renormalization group running our first scenario describes a 
crossover between two fixed points at 1/5 = and S — 0. The scale x corresponds 
to a dynamical infrared scale of the unknown fundamental theory. As it is lowered 
the flow of the couplings switches from the fixed point at = to the range 
of attraction of the conformal fixed point at 5* = 0. However, this latter fixed 
point is not stable with respect to the running of the gauge coupling which finally 
grows large and produces a new infrared scale m. This non-perturbative scale lifts 
the degeneracy in the fiat cosmon direction. This type of behavior is common in 
statistical physics when the crossover between two fixed points plays a role. 

Dark energy 

The cosmological dynamics of CQ is described in [19, 32] and best expressed 
after Weyl scaling with a constant reduced Planck mass = M^/Stt and a rescaled 
dimensionless scalar field $ = (p/M — 21n(x/m). In terms of the rescaled metric 
and cosmon field the effective action (20) reads [19] 

j d'x^ i^^[k\<f>)d'^<f>d,<j> -R\+ M^exp(-0) + M^F^^F/^-j (g) 

with ^ k'^{(/)) = S/4. The cosmon potential V = M^exp(— 0) vanishes for oo. 
In this language it is easy to see that the asymptotic behavior for increasing time 
corresponds to increasing and x- 

The cosmological evolution equations are best expressed in terms of the logarithm 
of the scale factor 

X = \na = — lia{l + z) , X = H. (9) 

The first equation expresses ^ the change of the energy density of homogeneous 
quintessence, ph, in terms of its equation of state Wh = Ph/ Ph, 

Ph + 3il + Wh)Hph = , ^^^ = -3{l + Wh). (10) 

^Eq. (4) corresponds to = 2E{(j)c — </>). Similar to eq. (7) this is a form of "leaping kinetic 
term quintessence" [19]. 

^ Wc have neglected the effects of the cosmon coupling to matter and radiation which are assumed 
to be very small in the present model. 
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Using Per = SM"^!!^ = Ph + Pm + Pr , dhipm/dx = —3 , dlnpr/dx = — 4 , Pr = 
{aeq/a)pm, eq. (10) can be translated into an evolution equation for flh = Ph/ Pen 
namely 

dQ,h 



dx 



"eq 



11) 



We recover the fixed points for w/j 1/3 if a = <^ a^q, and if » a^q. 

They are relevant for the cosmic attractors in the radiation and matter dominated 
epochs. 

As a second equation we need the time evolution of Wh for which we find 

dwh 



dx 



(1 - Wh) [k-\(t>W?,nn{l + Wh)- 3(1 + Wh)]. (12) 

The dynamics of the particular model enters through k~^{(f)) = 2/\/6. Our scaling 
of 4> is chosen such that at all times the value of 4> is directly connected with the 
value of the cosmon potential and therefore also with if^, 



(13) 



For given present values of the Hubble parameter Hq and the energy fractions in 
matter, radiation, and quintessence, f] 
expresses in terms of x, Q/j and Wh 



matter, radiation, and quintessence, ^Im , ^r'^ — deq^m , ^^^^ = 1 — — ^r^^ this 



4x-2 ln{Ho/M) - In 



aeq + e^n^^^ 



-In 



^1 - Wh)^h 



(14) 



The system of the two differential equations (11) (12) is therefore closed and can 
easily be solved numerically. 

For given E and Js = 0.05 we solve the flow equation (7). The "initial value" 
= ^{x = ^) is tuned such that at present = 0.73. We start for x{t = 0) = m 
with rather arbitrary values of Qh and w^. The late time behavior will not depend 
on this. In table 1 we present a few characteristic cosmological quantities, namely 
the present equation of state, w^^\ the age of the universe, t^^\ the fraction of 
dark energy at last scattering, the location of the third peak in the CMB 

unisotropies in angular momentum space, l^, as well as the normalization of the 
spectrum of density fluctuations, as, divided by the value which would obtain for 
a cosmological constant with the same fl^^^ = 0.73. Here the "present" values are 
characterized by the Hubble parameter reaching its present value with h — 0.71. 
We recall that smaller h and larger Q^^^ shift Is to larger values. 



E 


(0) 


^(0) /iQio^r 


"a 


/3 


/ (A) 


5 


-0.94 


1.31 


0.019 


779 


0.70 


12 


-0.99 


1.35 


0.0083 


794 


0.85 



^We use d(p/dx = fc-^ (617^(1 - V/ph)Y^^ and V/p^ = (1 - w/,)/2 
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Table 1: Characteristic cosmological quantities for two CQ models. 



We have chosen the cosmological parameters Q^^^ and h in accordance with the 
best fit of the KMAP-data [1] for the case of a cosmological constant (w^^^ = — 1). 
For E = 12 the values w^^^ = -0.99 and = 0.008 imply that for the CMB 
spectrum there is barely any difference as compared to the case of a cosmological 
constant with the same amount of dark energy [21]. We emphasize, however, that 
the value as/cTg^^ < 1 indicates that more power on large scales can be tolerated in 
order to achieve consistency with the galaxy surveys. This may make the KMAP 
data more compatible with a ^-independent spectral index. We have not tried to 
optimize the choice of and h. For = 5 a better agreement with the KMAP- 
data may be achieved by varying these parameters. We conclude that for large 
enough E our model is consistent with present cosmological observations ^. 

Cosmological phenomenology of crossover quintessence 

For a detailed comparison with cosmological observation one is interested in the de- 
pendence of the Hubble parameter on redshift, H{z), and related quantities. Rather 
than testing each particular model - as characterized here by the function k{(f)) - one 
wants to describe the function H{z) by a few parameters whose values may be fixed 
by future cosmological precision tests. We propose to concentrate for this purpose 
on the average equation of state (a; = — ln(l -|- z)) 



1 

Wh{x) = — / dx'whix'). (15) 

This function is linked very directly with cosmological observables and admits a 
simple parameterization. We believe that for the next generation of observations a 
quadratic formula will be sufficient (for z < 10^) 

Wh{x) = + {wl^"^ - + Ax{x - xis). (16) 

Xls 

Here the three parameters are the equation of state today, w^^^ , the averaged equa- 
tion of state at last scattering, w^f^^^ with xis — — In(llOO), and A. (In the future, 
one wants to measure the whole function Wh{x) as precisely as possible - it contains 
the complete information about the dynamics of dark energy.) 

Indeed, from Wh{x) one can easily construct the time history of flh according to 

ilj^ix) _ fif (1 + ae,) exp ( - 3.r»v,(.T;)) ^^^^ 



l-flh{x) 1 _ Q^O) 1 + Qeq exp(-x) 



^Adding to a term linear in 6 would further suppress O/, at early time and bring the phe- 

nonicnology even closer to the one for a cosmological constant (for given E). The prize to pay are 
unnaturally small values of 6i, similar to quintessence with an inverse power law potential [16]. 
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with tteq the scale factor at matter-radiation equahty. We can also obtain the equa- 
tion of state at any time 

Wh{x) = — {xwhix)). (18) 
The function H{z) is related to Wh{x) by 

^ = + .)3(i+-.(^)) + + ,)3 + a^^il + .)1 (19) 

where n^^^ = {1 - nf^^ ) / {1 + a^q). We note that wjj'*'' in eq. (16) is directly related 
to the fraction of dark energy at last scattering, via eq. (17). This quantity 

can be tested by the CMB-unisotropies [21]. For a given model of quintessence we 
determine A be equating the fit formula (16) with the model-computation of idh 
at X = Xis/2. For the models displayed in table 1 this yields reasonable fits for 
the range Xis < x < 0, with wf^'^^ = —0.221, A = —0.016 for E = 5 and wj^'*'' = 
—0.261, A = —0.013 for E = 12. We emphasize that the polynomial fit should not 
be used too far outside this range (e.g. for nucleosynthesis) since CQ typically results 
in a crossover from Wh near 1/3 for large negative x to a constant negative value of 
Wh for large positive x. However, for all observations of recent cosmology from the 
time of last scattering until now the proposed fit will do reasonably well. The fact 
that we choose Wh and a polynomial in ln(l + z) permits us to cover the whole range 
between last scattering and now. Supernovae observations at 2; ~ 1 can therefore 
be linked to structure formation and CMB within a simple parameterization. This 
would be impossible with a Taylor expansion of Wh{z) in z around z — 0. 

Cosmological history of fundamental couplings 

Under quite general circumstances quintessence or a time varying dark energy 
can be associated with a scalar field whose cosmological expectation value varies 
during the recent history of the universe [2]. Generically, this scalar field - the 
cosmon - may also couple to matter and radiation. As a consequence, the values 
of "fundamental constants" like the fine structure constant or the ratio between 
the nucleon mass and the Planck mass depend on the value of the cosmon field and 
therefore on cosmological time [9, 2, 17, 27, 28]. Recently, a time variation of the fine 
structure constant has been reported [12] from the observation of quasar absorption 
lines (QSO) at redshift z ^ 2. A typical value corresponds to 

^""'""^^ ^ = -7 • 10-^ (20) 

Similar observations infer ^ a substantially smaller value of lAoeml at lower redshift 
z < 0.7. From the Oklo natural nuclear reactor one obtains a typical bound [13], 

recent eollcetion of results on Aagm as well as the time variation of other fundamental 
constants can be found in [29]. 
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[14] 



|Aae^(^ = < lo-T (21) 



and analysis of the decay rate Re^^'^ — > Os^^"^ restricts [14] 



l^"^-^^ = °i^<3.10- (22) 



As an obvious question one may ask if an increase of |Aaem| by almost two orders of 
magnitude between z = 0.13 and z = 2 can be reasonably explained by quintessence 
models. In this note we demonstrate that this can indeed be the case in a class of 
models of "crossover quintessence" proposed recently [11]. This contrasts to constant 
rates dAaem/dz = const, or dAaem/dt = const, for which the values and bounds 
(20)-(22) clearly are in discrepancy. The deviation from constant rates also affects 
strongly the comparison of eq. (20) with bounds from nuclear synthesis [30], [31] and 
CBM [31]. 

We are aware that the bound from the Oklo natural reactor is subject to sub- 
stantial uncertainties due to possible cancellations between effects from the variation 
of aem and other "fundamental constants" like the mass ratio between pion and nu- 
cleon mass m,T^/m,n or the weak interaction rates. Also the QSO observations need 
further verification and investigation of systematics. Nevertheless, it is well con- 
ceivable that we can get a reasonable "measurement" of the function Aaemi^) ^ 
not too distant future. Furthermore, the derivative dA.aem{,z) / dz at z = can be 
related [32] to precision tests of the equivalence principle [15]. The aim of this 
note is to present a sample computation how the knowledge of Aaem{z) can be used 
to constrain models of quintessence. For this purpose we will use the value (20) 
as a "benchmark value" which fixes the strength of the cosmon coupling to matter 
and radiation. The time history Aaem{z) can then be related to the time history of 
quintessence. 

In particular, we may investigate the ratio 

R = Aaem{z = 0.13)/AQ;em(-2 = 2) (23) 

A bound R < 1/70 strongly favors [32] quintessence with a time varying equation 
of state Wh = Ph/ Phi where the value of (1 -I- Wh) at present is substantially smaller 
than for 2; = 2. A logarithmic dependence of Aaem on the cosmon field x leads to 
a bound for the present equation of state [32] 

< -0.9. (24) 

In this note we investigate how the detailed dependence of on x influences the 
time history of Aaem{z). We argue that reasonable /9-functions for the "running" 

^°Sec [27], [28] for a different relation which gives a much smaller violation of the equivalence 
principle. 
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of a grand unified coupling with Inx may lead to nonlinearities in Q;em(lnx). Those 
are quantitatively important, influencing R within a factor two. The overall picture 
remains rather solid, however: a small value of \R\ requires a quintessence model 
where the evolution of the cosmon field has considerably slowed down in the recent 
history as compared to a redshift z ^ 2. This feature is characteristic for crossover 
quintessence and will not be shared by arbitrary quintessence models. In particular, 
a constant equation of state with Wh independent of z over a range < z < 3 will 
have severe difficulties to explain R < 0.1. This demonstrates how a measurement of 
Aaem{z) could become an important ingredient for the determination of the nature 
of dark energy. 

Running couplings 

In a grand unified theory the value of the electromagnetic fine structure constant 
depends on the gauge coupling at the unification scale, ax, the ratio between the 
weak scale and the unification scale, M^/Mx, and particular particle mass ratios 
like rrin/Mw or rrie/Mw for the nucleons and electrons. Since the nucleon mass is 
determined by the running of the strong gauge coupling similar dependencies arise 
for the ratio m„/Mx. Combining the relevant one loop formula and assuming for 
simplicity x-independent ratios rrib/Mw etc. for the heavy quarks and m^/me etc. 
for the leptons one finds [32] 

1 1 22 



aera{z) Q!em(0) 7 \ax{z) ax(0) 



2l7r (Miy/m„)(0) tt (Mw./me)(0) 



and 



(mjM)(0) 11 \a,„,{z) ) ^ ^ 

^ ^ {Mw/mn){z) ^ 2_ {Mw/m,){z) ^ {Mx/M) {z) 



33 {Mw/rrinm 11 {Mw/m,){0) (Mx/M)(0)' 

We will make here the simplifying assumptions that the effects due to the varia- 
tion of Mw/mn and Mw/rUe are subleading, as suggested by the more complete 
investigations in [32] (in the linear approximation). 

In the approximation of time- independent Mw/rrin and Mw/rrie the dependence 
of the fine structure constant on redshift can be related directly to the /3-function 
for the grand unified gauge coupling (2). For /^aem{z) — aem{z) — Q;em(0) one finds 
(x = - \n{l + z)) 

Aaem{z) 22a 

e 



'-Aaxix) (27) 



^We have eliminated the dependence of Ofem on Mx/M^y in favor of Mty/rn„. 
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where ax{x) obtains as a solution of 



dax ^ <91nx PaioxJ) / , r , . 

^ = /5«^ = ^^v/^iMrT^. (28) 

One observes that the ^-dependence of Acte^ is particularly weak in the region of 
large S and Wh close to —1. This will be our explanation why R = Aaemi^ — 
0.13)/Aaem{z = 2) is considerably smaller than expected from a simple extrapo- 
lation linear in z. Crossover quintessence is precisely characterized by large 6 and 

~ — 1 in the region of small z, whereas \/l + Wh/ \/S is considerably larger for 
intermediate 2; 1.5 — 3. One also observes that for constant (Jl/j(l -|- Wh)/SY^'^ 
and r] — (5a/ oc the dependence of A.aem{,z) is logarithmic in 1 + 2;. This is crucial 
for large z as for nucleosynthesis, where Aaem turns always out to be much smaller 
than expected from a linear extrapolation in redshift or time. 

A quantitative estimate needs information about Pa - the form of the latter 
may induce further nonlinearities. Prom our association of the scale m with the 
scale where ax grows large the /3-function (2) would lead to a decrease of ax for 
decreasing z if the last term ~ S{5) is neglected. This would result in a positive 
/\aem{z = 2), in contrast to the negative QSO value (20). Starting with large 
ax{x = ^) (we take axix = ^) = 100 for definiteness) one would actually find 
values of \A.aem{z — 2)| much smaller than the QSO value if 62 exceeds 0.1. This 
demonstrates that a fixed point behavior for the running of ax (and similar for other 
couplings) provides a powerful explanation why the time variation of "fundamental 
constants" is actually a small effect! 

There is, however, no reason why the fixed point value a^, should be precisley 
independent of 5. We parameterize a possible 5-dependence by adding the last term 
in eq. (2) with 

This leads to a crossover of the effective fixed point a^{5) between two values for 
small and large 5. For 65 > , > the value of a* (5) increases with increasing 5 
and this will induce /S.aem{z = 2) < 0. At this stage the form of (5a is parameterized 
by four parameters 62, ^'4, -Be, Since the influence of the term S{5) on the location 
of the fixed point is small the ratio 62/^*4 is fixed by eq. (3). 

Time variation of the fine structure constant 



At this point our model is fixed and we can now compute the detailed cosmolog- 
ical evolution of the fine structure constant. The overall normalization of the time 
variation strongly depends on h^. As a benchmark we use eq. (20). Therefore, for 
given 62 and we fix feg such that A.aem{z = 2) = — 7 - 10~^. The shape of A.aem{z) 

^^We choose S{5 = 0) > such that for 65 > no new fixed point for ax is produced in the 
region of large ax- The value S{5 = 0) = 6 is unimportant. 
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depends now on 62 and J a as well as on on In table 2 we show the values of 
Actem at z = 0.13, 0.45, 1100 and 10^° for several values of (62, Ja, E). 



model 


z = 0.13 


z = 0.45 


z = 1100 


z = lOio 


R 


V 


A 


-6 ■ 10-8 


-2.6 ■ 10-^ 


-6.5 ■ 10-5 


-7.8 ■ 10-5 


0.0085 


5.1-10-14 


B 


-6.2 ■ 10-8 


-2.7-10"^ 


1.3 ■ 10-4 


3.1 ■ 10-'^ 


0.0088 


4.9-10-14 


C 


-7.8 ■ 10-8 


-3.5 ■ 10-^ 


-2.5 ■ 10-5 


-2.8 ■ 10-5 


0.011 


8.9 - 10-14 



Table 2a: Time variation of the fine structure constant Aaem(-2)/ctem for various 
redshifts z, for E = 12. The models for j3a are (A) : = 6, 62 = 0.2, fog = 
0.84, (B) : J„ = 6,62 = 0.05,^6 = 1.12, (C) : J„ = 1,62 = 0.2,66 = 0.15. For all 
models Aaem(-s = 2)/Q;em = ~7 - 10"^. We also show the differential acceleration rj. 



model 


z = 0.13 


z = 0.45 


z = 1100 


z = lOlo 


R 


V 


A 


-1.4 - 10-^ 


-7.1 - 10-^ 


-4.4 - 10-5 


-5.7-10-5 


0.020 


4.3-10-14 


B 


-2 - 10-^ 


-9.3 - 10-^ 


1.7- 10-4 


3.6- 10-3 


0.028 


6.2 - 10-14 


C 


-2 - 10-^ 


-9.9 - 10-^ 


-2.1 - 10-5 


-2.3 - 10-5 


0.029 


8.7-10-14 



Table 2b: Same as table 2a, for E — h. The models for (3a are (A): Jq, = 6, 62 = 
0.2, be = 0.385, (B): J„ = 6, 62 = 0.05, = 0.783, (C): = 1, &2 = 0.2, feg = 0.068. 



Consistency with the Oklo bound (21) and the Re-decay bound (22) is found 
for large in particular for large J a- It becomes clear that a precise knowledge of 
^aem{z) could be used as a probe for the dynamics of quintessence. For example, 
it seems very difficult to achieve R = \Aam{z = 0.13)/AQ;em(2; = 2)| < 1/70 
for small values of E, quite independently of the precise form of j3a- The present 
investigation confirms that a small value of R requires (l + wf^) -C 1. This property 
of the equation of state seems to be favored by the cosmological tests (cf. table 1). 
We conclude that the bound (24) can be justified by our nonlinear analysis for CQ. 
The generic result that w^^^ should be close to —1 does not depend on the precise 
values of 62 and J^. Increasing 62 beyond 62 = 0.2 (say 62 = 1) induces almost no 
change. On the other hand, for 62 substantially smaller than 0.05 we come to a 
region where the fixed point ax,* has not yet been reached with high precision and 
Auemi^ — 2) is positive. For 62 = 0.05 we have the interesting effect that Aagmi^) 
reaches a minimum near z = 2, changing sign at somewhat higher redshift. Only for 
the large value E = 12 the functional dependence of Aaem{z) deviates substantially 
from a linear behavior in z in the range 0.5 < 2; < 2, as seems to be required by 
observation. For the models (A) (B) (cf. table 2) one finds a strong "jump" of Actem 
by a factor of 5 between z — 1 and z — 2. 

Let us next turn to the high values z ~ 10^^ characteristic for nucleosynthesis 
and ask if the values for Aagm quoted in table 2 are consistent with observation. 

^^In our case an increase of |/3a| for smaU z partiafly cancels the decrease of (f2/j(l + Wh)/Sy/^. 
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For nucleosynthesis, a change in nin/ M modifies the clock, i.e. the rate of decrease 
of temperature in time units ~ Furthermore, a change in ctem affects the 

proton to neutron mass difference. Let us assume that the time- variation of m„/ M 
is dominated by (3^ via the change of the strong gauge coupling (and therefore the 
ratio Aqcd/Mx)- Including only the first term in eq. (26) yields 

A{mn/M) _ TT Aaem ^ gg ^ AOem ^g^^ 

A decrease of m„ relative to M is equivalent to an increase of M at fixed m„ (and 
thereby to a decrease of Newton's constant). In turn, this decreases the Hubble 
parameter for a fixed temperature T (measured in units of m„). The net effect is 
equivalent to the substraction of a (fractional) number of neutrino species. We can 
thereby turn the limits on the effective number of neutrino species N,^ into a limit 
on Aaem/(^em Sit z ^ 10^°. Assuming that the primordial helium abundance Yp is 
in agreement with the value for time independent couplings within an uncertainty 
\AYp/Yp\ < 8 • 10"^ and using AlnYp = |Aln(m„/M), yields the bounds 



and 



|Aln(mn/M)| < 0.025 (31) 



I ^1 < 6.4 -10"^ (32) 



We emphasize that an effective neutrino number < 3 would be a clear signal for 
a time variation of couphngs! We also observe that bounds on m„/M as (31) can 
be interpreted equivalently as bounds on the time variation of Newton's constant 
(with m„ fixed). 

Wc may compare the effect of the change of m„/M with the effect of the change 
in the proton-neutron mass difference dm = m„ — Trip which contributes AlnYp ^ 
iAln(m„/M) — A In (5m. Keeping the ratios of up-and down quark masses to m„ 
fixed one has [34] A In 5m ~ — 0.6A In ctem- The contribution to Alnl^ is far less 
than the contribution from Aln(m„/M). This explains why the bound (32) is more 
restrictive than several previous bounds [30], [31] concentrating on the effect from 
A In 5m. We are aware that the variation of Aln(m„/M) may be affected by the 
neglected terms in eq. (26). If there are no particular cancellations, however, the 
bound (32) should remain of the same order of magnitude. 

The size of Aaem at nucleosynthesis and even the sign depend critically on 62 
(cf. table 2). If we adopt the bound 

] ^1 < 10"^ (33) 



-'^''This relation differs from [33] where Mw/Mx is kept fixed instead of fixed irin/Mw in the 
present work. The quantitative difference is not very important, however. 
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we infer that 62 should be larger than 0.05. It is striking that the values of Aaem/ctem 
at nucleosynthesis are very far from an extrapolation with din aem/dt — const, or 
even from d In ccem/ dlna — const.. The nonlinearity in d In ccem/ 9 In a is mainly due 
to the attraction of ax to a 5-dependcnt fixed point. Both for E = 5,12 the models 
with sufficiently large 62 meet the bounds (32) or (33) (model A). They are also 
consistent with bounds from the CMB unisotropies [31]. 

Precision experiments with atomic clocks measure the variation Aaem/cxem per 
year 

-[yr \ = -H[yr 



ocem d\n(l + z) 

Here the Hubble parameter is expressed in units of yr~^ and we use the observation 
that for our models dlnaem/dx is essentially constant in the range < 2; < 0.13. 
Measuring the values quoted in table 2 would need a considerable improvement of 
the present accuracy [35] dem/ctem = (4.2 it 6.9) ■ lO^^^yr^^. Using eq. (30) we 
obtain a similar formula for the variation of Newton's constant in the range of low 
redshift, 

^ = 78.2^ = -4.25 • io-8 Ac.em(^-0.13) ^ ^3^^ 

safely below the present bounds. 

In our approximation where the dominant field dependence of the various cou- 
plings arises from (3^ we can also estimate the size of the differential acceleration 77 
between two test bodies with equal mass but different composition. One finds [32] 

rj ^5 -10-^ (J^y S-'ARz (36) 

where for typical experimental tests of the equivalence principle ARz — AZ/{Z + 
N) 0.1. Here Pa/ctx and 6 have to be evaluated at z = 0. Results for rj are also 
displayed in table 2. We conclude that our CQ-model is compatible with the present 
experimental bounds [15] < 3 ■ 10^^^ for most choices of Pa- Nevertheless, the 
models with Jq, = 1 show that already now the tests of the equivalence principle can 
help to discriminate between different models! It is obvious that further improve- 
ments of the accuracy of tests of the equivalence principle should either confirm or 
reject the interpretation of the QSO-results within our CQ-model. 

In summary, we have found models of crossover quintessence that can explain the 
QSO value for the time dependence of the fine structure constant and are neverthe- 
less compatible with all observational bounds on the time variation of couplings and 
tests of the equivalence principle. This explanation does not involve any particular 
cancellation of effects from the variation of different couplings. The parameter set 
E = 12, 62 = 0.2, Jq, = 6 may serve as an illustration. The same models do also very 
well with all present cosmological tests! 
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Variation of electron-proton mass ratio 



Let us finally ask if our assumption that Accem is dominated by the change of ax 
can be justified. In addition to the arguments presented in [32] we may also have 
a look at the recently reported result [36] on the variation of nie/mp in the range 
^ = 2.3-3 

The deviation from zero is statistically not very significant. Nevertheless, it is inter- 
esting to know if such an effect would influence the time history of the fine structure 
constant and the tests of the equivalence principle. We therefore investigate what 
would be the effect on Actem of a variation A ln(me/m„)(z = 2) = — 6 - 10~^. Keeping 
Mw/rUn fixed would yield a contribution (cf. eq. (25)) 



A In — a; -2.8 • 10~\ (38) 

Otem TT J 

This amounts only to 4% of the QSO value (20). Similarly, the contribution (cf. eq. 
(26)) to 

m i^-f;;--^; . - Aa m (^) . 1.1 . 10- (39) 

{mjM){z = Q) 11 \mn) ^ ' 

is only 4% of the effect due to Aagm, i-e. 

,^(m„/A?)(. = 2)^ . ^"°"(^ = ^) ^-2.7.10-". (40) 

Neglecting those contributions seems therefore well justified. Comparison of eqs. 
(40) and (37) incidentally shows that the value (37) corresponds roughly to the size 
of the effect that would be expected in absence of particular cancellations - with 
ln(me/M)(^ = 2)/ln(me/M)(;2 = 0) -3.3 • 10"^. The effect of Aln(me/m„) on 
the tests of the equivalence principle is more sizeable. Adding this effect to the 
differential acceleration multiphes r] in eq. (36) by a factor l^Q [32] 

~ d(me/mn) A(me/m„) , , 

Q^IO ^ -10 V ^/ '" {z^2)^6, (41) 

daem [z=0 ^Olem 

where the quantitative estimate neglects a possible 2;-dependence of Q. Multiplica- 
tion of the ?7- values in table 2 by a factor 7 shows that for part of the models a value 
A ln(me/mp)(z = 2) = —6 ■ 10"^ enters in conflict with the tests of the equivalence 
principle. The models with E = 12 and Jq, = 6 may be considered as borderline 
in view of the possibility of partial cancellations with effects from the variation of 
other mass ratios. 
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The future universe 



If we live today a crossover period of our universe, to which future will the 
transition lead? In our model, the answer depends in important aspects on the 
behavior of (3s for large 5. The scenario that we have discussed first is characterized 
by ever increasing 6 (or the approach to a fixed point (1/5)* = 0). In this case 
Q/j will grow towards one and Wh approaches —1 in future times. In these respects 
the universe will resemble asymptotically a universe characterized by a cosmological 
constant. Nevertheless, if k{(j)) remains finite for all finite 0, the dark energy density 
Ph will not approach a constant since W/i + 1 will always remain positive, approaching 
zero only for (p ^ oo. An asymptotic behavior (3^ = D6 for large 6 corresponds 
to k{(f)) — /c^exp 1^(0 — 0c)}- This is the asymptotic behavior of inverse power 
law quintessence [16]: rescaling the cosmon kinetic term leads to a potential 
V{(p) ~ <^~" , a = 4/D. In leading order one finds l + w^ = 1/(3A;^). The potential 
energy decreases logarithmically with a, V ~ (lna)~^/^ and we infer the asymptotic 
behavior ~ (lna)^/^/a"^, ilh = ^ — ^m- 

As a second scenario we discuss the alternative that (3s has a zero at some finite 
value 5*. As a natural choice we take — 6, noting that for 5 = 6 the kinetic term 
for X in 6q. (1) vanishes, and consider 

Ps = Ed' (l - • (42) 

Our cosmological model describes now a crossover between the two fixed points 
— and — 6, being relevant for the early and late universe, respectively. The 
solution 6{x) to eq. (42) obeys 

X 1 

(43) 



l + |ln(i-|) Elnixc/X) 



and approaches the solutions of eqs. (4) or (7) for small S. The asymptotic cosmology 
for large t is different, however. It corresponds to "power law inflation" where the 
equation of state approaches the asymptotic value — —7/9, with 

36 

5 = 6-— exp[-3£;(0 - (f)c)]. (44) 

The Hubble parameter decreases ~ t~^, with a ~ t^. Comparing ~ a~'^^,Pm ^ 
we conclude that Q/j — > 1, Qrn ~ a~^/^. As for the first scenario, the future 
universe will be completely dominated by dark energy. However, the dynamics of 
the dominant quintessence component differs between the two scenarios. 

^^For our ansatz (7) with D = E/Js we find very small a, i.e. for E = 5,Js = 0.05 follows 
a = 0.04. 
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We have displayed the characteristic cosmological observables in table 3. In 
order to demonstrate the spread in this type of models we have chosen h — 0.65 and 
= 0.7 and also display smaller values of £■ as compared to table 1. 



E 


(0) 


tW/10i"?/r 




k 


/ (A) 


2 


-0.65 


1.31 


0.045 


768 


0.46 


5 


-0.81 


1.37 


0.019 


795 


0.69 



Table 3: Characteristic cosmological quantities for CQ-models with fixed point at 
5* = 6. 

We also show in table 4 the time variation of aem, using in eqs. (2) (29) — 
and tuning ftg such that Aa^rniz = 2)/Q!em = ~7 • 10"^. 



E 


b2 


h 


Aaem{z = 0.13) 


2 


0.2 


0.011 


-1.1 ■ w-^ 


5 


0.06 


0.022 


-1.8 ■ 10-6 


5 


0.2 


0.01 


-9.5 • 10"^ 



Table 4: Variation of Aaemiz — 0.13)/Q;eTO for fixed Aaemi^ — '^)/'^em = — 7 • 10 ^. 
The models for use = 0. 

It becomes obvious that the bound (21) is not obeyed. Consistency with the Oklo- 
natural reactor results would therefore require a substantial cancellation between 
effects from the time variation of different couplings. 

Testing the model 

Within a rather wide class of quintessence models we have seen how a com- 
bination of the QSO-result (20) with the bounds on «em at lower redshift (21) 
(22) restricts the allowed parameter space. The present status of observations 
favors a rather sharp crossover for quintessence, where the equation of state has 
shifted from a moderately negative value at z 2 to a value very close to — 1 
for the recent cosmological epoch. An example is given by the effective action (1) 
with running couplings determined by the /5-functions (7) (2) (29), and parameters 
E = 12, Js = 0.05,62 = 0.2,64 = 8,66 = 0.84, J« = 6 (first line in table 1, model 
A in table 2a). Without the need of excessive tuning of parameters or invoking 
cancellation effects from the time variation of different couplings, this demonstrates 
the existence of a field theoretical model that is consistent with all data on the 
cosmological evolution and the time variation of couplings. 

What will be the tests of this model? First of all, the QSO result should be 
confirmed by new independent data. Hopefully this will allow for more precise re- 
strictions on the shape of the function A In Q;em(^) which can be directly compared 
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with the prediction of this model. Second, an improvement of the accuracy of the 
tests of the equivalence principle by an order of magnitude should lead to a di- 
rect detection of the new interaction mediated by the cosmon field. In view of its 
important role in present cosmology we may call this fifth force (besides gravity 
electromagnetism, weak and strong interactions) "cosmo-dynamics" . Third, cosmo- 
logical tests should distinguish crossover quintessence from a cosmological constant. 
The main signature in this respect is probably not the equation of state of dark 
energy in the most recent epoch (say for < z < 0.5). The value of wf^ in table 1 
shows that this quantity may be quite near to the value for a cosmological constant 
(i.e. wf'^ = — 1). If the QSO-results on the time variation of the fine structure 
constant are to be explained by quintessence this requires a substantial evolution 
of the cosmon field in the epoch near 2; = 2 (as compared to the rate of evolution 
today). We therefore expect that Wh{z) increases with increasing z, leading to ^ui^) 
substantially larger than for a cosmological constant during structure formation and 
last scattering. Test of structure formation and precision tests on the CMB may be 
able to measure the value of as/cr^^ and Vl^^^^ in table 1. Finally, an improvement 
of the accuracy of atomic clocks by three orders of magnitude would replace the 
bounds (21) (22) within an experimentally well controlled setting. 
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